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Abstract 

We report on measurements of low frequency noise in a single electron transistor from a few Hz 
up to 10 MHz. Measurements were done for different bias and gate voltage, which allows us to 
separate noise contributions from different noise sources. We find a 1/f noise spectrum with two 
Lorentzians superimposed. The cut-off frequency of one of the Lorentzians varies systematically 
with the potential of the SET island. Our data is consistent with two single-charge fluctuators 
situated close to the tunnel barrier. We suggest that these are due to random charging of aluminum 
grains, each acting as a single electron box with tunnel coupling to one of the leads and capacitively 
coupled to the SET island. We are able to fit the data to our model and extract parameters for 
the fluctuators. 

PACS numbers: 07.20.Mc, 73.23.Hk 
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I. INTRODUCTION 



Single electron transistors (SETs)^ - 1 ^ are extremely sensitive electrometers, with demon- 
strated charge sensitivities of the order of /ie/v / Hz 1,5 . Due to their high charge sensitivities 
they have found a large number of applications in research, for example, SETs are used to 
detect nano-mechanical oscillators^, to count electrons^" 22 and to read out qubits&i^i 

The fundamental limitation for the sensitivity of the SET is set by shot noise generated 
when electrons tunnel across the tunnel barriers^- 5 . Shot noise was observed in a two 
junction structure (without gate) by Birk et al. 3 . 

However, there are two other types of noise which are limiting the charge sensitivity in 
experiments. At high frequencies, above 1 MHz, the sensitivity is limited by the amplifier 
noise. At low frequencies the sensitivity is limited by 1/f noise which is due to background 
charge fluctuators near the SET. A collection of several ffuctuators with different frequencies 
leads to a 1/f spectrum. In several cases it has been observed that there is a background 
of 1/f noise and a single or a few more strongly coupled fluctuators, resulting in random 
telegraph noise, which in the frequency spectrum leads to Lorentzians superimposed on the 
1/f background 2 ^. 

Understanding the nature of the 1/f noise is also very important for solid state qubits 
since 1/f noise strongly limits the decoherence time for these qubits. It has been suggested 
by Ithier et al. that the charge noise has a cut-off at the frequency of the order of 0.5 MHz^. 

Even though there have been many efforts^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^'^ 1 ^ to reveal the physical 
origin of the background charge fluctuators the nature of these fluctuators is still unknown. 
It is not even clear where these fluctuators are located. The charge fluctuators can be located 
either in the tunnel barrier or outside the barrier but in close proximity of the junction. 

The role of the substrate has been examined in several experiments^ 1 ^ 120 -. However, those 
experiments did not show a strong dependence of the noise on the substrate material. The 
barrier dielectric has been proposed as location of the charge traps by several groups^ 2 ^ - 1 ^. 

Several groups have shown that the low frequency noise at the output of the SET varies 
with the current gain (dl/dQg) of the SET and that the maximum noise is found at the 
bias point with maximum gain 2 ^ 1 ^^ -. This indicates that the noise sources acts at the input 
of the device, i.e. as an external fluctuating charge. A detailed comparison of the noise to 
the gain was done by Starmark et a/.—. All above mentioned experiments were performed 
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with conventional SETs by measuring current or voltage noise at relatively low frequencies. 
i.e. below a few kHz. 

In this work we have measured low frequency noise in the Single Electron Transistor 
which has demonstrated a very high charge sensitivity 5 , by using the Radio Frequency Single 
Electron Transistor (rf-SET) technique 25,3 *. This allowed us to measure low frequency noise 
of the reflected voltage from the rf-SET in the range of a few hertz up to tens of MHz, 
and due to high charge sensitivity we were not limited by the amplifier noise. We find two 
Lorentzians superimposed on a 1/f spectrum and that the noise in the range50 kHz-1 MHz is 
quite different for positive and for negative gain of the transistor. By analyzing the bias and 
gate dependence of the noise we argue that the noise in this frequency range is dominated 
by electron tunneling to an aluminum grain, which acts as a single electron box capacitively 
connected to the SET island. 

The paper is organized as follows. In section [III we describe a model for the low frequency 
noise, which allows us to separate contributions from the different noise sources. In section 
IIHI we describe the experimental details of our measurements. Section IIVI is the main part 
of this paper and contains the experimental results. Finally in section |V] we describe our 
model for the nature of the low frequency noise in our SETs. 

II. LOW-FREQUENCY NOISE MODEL 

We start by analyzing the different contributions of the measured noise. What we actu- 
ally measure is the reflected voltage from the tank circuit in which the SET is embedded. 
The rf-SET tank circuit is a series LC circuit with inductance L and capacitance C. The 
power spectral density of the reflected voltage can be decomposed into the following terms 
originating from charge noise, resistance noise, shot noise and amplifier noise according to 
the following equation— 1 ^: 

Sm = ( Sqs >^ + ("J^") Sri ^ + (a^) Sr2 ^ + Sshat + SAm P i -> W 

where Q g is the charge at the input gate and R\ 2 are the tunnel resistances of the two 
junctions. Here we neglected the higher order terms and possible correlation terms between 
the charge noise and the resistance noise. In the case when the charge fluctuator is located 
in the tunnel barrier the correlation term may not be negligible. By measuring the noise for 
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different bias points having different gains (i.e. different d\v r \/dQ g ) it is possible to extract 
information on whether the noise is associated with charge or resistance ffuctuators. 

We have designed the matching circuit for the rf-SET to work in the over-coupled regime, 
in order to have a monotonic dependence of the reflection coefficient as a function of the 
SET differential conductance. This regime corresponds to the matching condition when 
the internal quality factor of the rf-SET tank circuit (Qi n t = y/L/C/R) is larger than the 
external quality factor (Q ext = Z / y/L/C), where Z = 50 f2 is the characteristic impedance 
of the transmission line, and R is the SET differential resistance. 

We have theoretically analyzed the reflected voltage from the rf-SET as a function bias 
and gate voltages applied to the SET. The reflected voltage characteristic of the rf-SET 
can be calculated by the orthodox theory^ using a master equation approach. From this 
theory we can also calculate the derivatives in eq. (JT]) of the reflected voltage with respect to 
variations in gate charge and in the resistances of the SET tunnel junctions. 

Figure Ufa) shows the sensitivity of the rf-SET to charge fluctuations as a function of the 
bias and gate voltages. The charge sensitivity is a symmetric function around the SET open 
state (Q g = 0.5 e), and has maxima close to the onset of the open state. 

The sensitivity of the SET to resistance fluctuations in the first tunnel barrier is shown in 
the figure DJb). The sensitivity to resistance fluctuations in the second barrier (not shown) 
is identical to figure [T](b), except that it is mirrored along the SET open state (Q g = 0.5 e). 

By operating at different bias and gate voltage we can choose operation points where the 
noise contribution from the different derivatives dominates, and it is possible to distinguish 
charge noise from resistance noise. 

III. EXPERIMENT DETAILS 

The samples were fabricated on oxidized silicon substrates using electron beam lithogra- 
phy and a standard double-angle evaporation technique. The asymptotic resistance of the 
measured single electron transistor was 25 kQ. The charging energy Eq = e 2 /2Cs ~ 18±2 K 
was extracted from the measurements of the SET stability diagram of the reflected signal 
with frequency / = 350 MHz versus bias voltage and gate voltages. From the asymmetry 
of the SET stability diagram we could also deduce that the asymmetry in the junction 
capacitances was 30%. 
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FIG. 1: (color online) Calculated derivatives from equation ([T]) as a function of bias voltage V, 
and gate charge Q g = C g V g . The derivative of the reflected voltage from the rf-SET with re- 
spect to (a) the charge fluctuations (d\v r \ / dQ g ) 2 ; (b) the resistance fluctuations in the 1st junction 
(d\v r \/dRi) 2 . Sensitivity to the resistance fluctuations in the second junction has a mirror sym- 
metry, along SET open state (Q g = 0.5 e), with the sensitivity to the resistance fluctuations in the 
first barrier. 

The SET was embedded in a resonant circuit and operated in the radio frequency 
modo^^i. The bandwidth of the setup was approximately 10 MHz limited by the qual- 
ity factor of the resonance circuit. The radio frequency signal was launched toward the 
resonance circuit and the reflected signal was amplified by two cold amplifiers, and then 
downconverted using homodyne mixing. The output signal from the mixer containing the 
noise information was then measured by a spectrum analyzer. The sample was attached to 
the mixing chamber of a dilution refrigerator which was cooled to a temperature of approxi- 
mately 25 mK. All measurements were performed in the normal (nonsuperconducting) state 
at a magnetic field of 1.5 T. 
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FIG. 2: (color online) (a) The SET stability diagram (dI/dQ g ) measured at a temperature of 
25 mK. Horizontal black line corresponds to the bias voltage, where transfer function I(Q g ) (see 
panel (b)) was obtained. The points A/ ((o) for negative and (+) for positive bias) where 1 = 
inside the Coulomb blockade region. The points B and D (*/•) correspond to maximum posi- 
tive/negative gain, dl/dQg that is where d 2 I /dQ 2 g = 0. The measurement points C, close to the 
current maximum current, correspond to 8I/dQ g = 0, marked with (+) at negative bias and with 
(o) at positive bias, (b) SET transfer function I(Q g ) measured for the bias voltage V = — 0.4mV. 
In the stability diagram this measurements is shown as a black line. 

We have performed the noise measurements for different gate voltages and different bias 
voltages. Due to experimental problems we have performed measurements mostly for neg- 
ative biases. The sample shows very high charge sensitivity of the order of 1 /ie/vTIz. For 
a detailed description of the sensitivity with respect to different parameters (see ref.~). A 
small sinusoidal charge signal of 7.3 • 10~ 4 e rms with a frequency of 133 Hz was applied to the 
gate, which allowed us to calibrate the charge sensitivity referred to the input of the SET. 

Before each measured spectrum, we have employed a charge-locked loop^, with a help of 
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a lock-in amplifier. The first (dI/dQ g ) or the second (d 2 I /dQ 2 g ) derivatives of the current 
were used as an error signals for stabilization of the gate point to compensate for the slow 
drift at the current maximum points (C in fig|2^b)), or at the maximum gain points (B and 
D in figfS^b)) respectively. Each noise spectrum is however measured with the feed back 
loop turned off. 

IV. EXPERIMENTAL RESULTS 

In order to separate the contributions from different noise sources we have performed 
measurements at specific points. The measurement points are shown in fig. [2j At point A 
[Q g k Oe) there is an almost a complete Coulomb blockade with a zero current and zero 
gain (d\v r \/dQ g = 0). Here the derivatives of the reflected voltage with respect to resistance 
fluctuations in the tunnel barriers (d\v r \/ dR\p = 0) are also zero (see figfjjb)). At this 
point we see only amplifier noise — different curves for different bias voltages show the same 
noise level convincing us that we see only amplifier noise. These measurements thus serve 
to calibrate the noise of the amplifier. 

The measurements at points B and D, correspond to the requirement of maximum current 
gain (max \dI/dQ g \) and therefore also high \d\v r \/dQ g \, diagonals in figur^a). At these 
points there are contributions from all the noise sources (see fig. [T]), but since the absolute 
gain and the current are very similar at the points B and D, these measurements can be 
compared. 

We have also measured noise at the points C (Q g ~ 0.5 e) close to the maximum of the 
current transfer function (\I(Q g )\). At this point the gain of the reflected signal (d\v r \/dQ g ) 
is quite low, but the shot noise could be high. The derivatives of the reflected voltage with 
respect to resistance fluctuations in the tunnel barriers (d\v r \/dRi^) are small but not equal 
to zero (see fig. H^b)). 

The noise of the reflected voltage at the fixed bias point and for the different gate points 
are shown in fig. [3](a). We start by subtracting the amplifier noise and then we compare 
the noise measured at the different points. Comparing the noise measured at points B and 
D, where the current gain has a maximum, with the noise measured at point C close to the 
maximum of the transfer function we see that the noise at the point C is substantially lower, 
even though the current is higher. From this we draw the conclusion that the difference is 
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FIG. 3: (color online) (a) Power spectral density (PSD) of the reflected voltage measured in the 
points A (black curve); B (blue curve); C (green curve); D (red curve), (b) Normalized noise in 
the points B (blue curve); D (red curve). The black continuous lines are a fits to the measured 
PSD with a sum of two Lorentzian functions. 

not due to the shot noise. 

Comparing the noise spectra measured at the points B and D it is necessary to note that 
both the currents and the gains are very similar at these points. Both spectra B and D have 
a 1/f dependence at low frequencies with two Lorentzian shoulders at higher frequencies. 
At frequencies above 30 kHz the noise at point D drops well bellow the noise at point B. At 
300 kHz the difference is a factor of 5. 

We have fitted the noise spectra at points B and D to a sum of two Lorentzian functions. 
The results of these fits are shown in figure [3](b). From these fits we can extract the cut-off 
frequency and the level for each of the two Lorentzians. 

The low- frequency Lorentzian has a cut-off frequency of the order of 1 kHz. The cut-off 
frequency is the same for both slopes (dI/dQ g <0), and it does not show any bias or gate 
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dependence within the accuracy of our measurements. 

In contrast, the high-frequency Lorentzian has a cut-off frequency above f co > 50 kHz, 
with a strong dependence on the bias and gate voltages. The bias dependence of the 
Lorentzian cut-off frequency for the positive (dI/dQ g > 0) and negative (dI/dQ g < 0) 
slopes are shown in figure S](a). 

For the negative slope (dI/dQ g < 0)(blue points in figure H]( a)), the cut-off frequency 
remains practically constant (f co ~ 50 kHz) for negative biases. Close to the zero bias voltage 
the Lorentzian cut-off frequency switches to a higher frequency (f co ~ 150 kHz). For the 




FIG. 4: (color online) (a) The bias dependence of the cut-off frequency for the high frequency 
Lorentzian. (b) The gate dependence of the cut-off frequency for the high frequency Lorentzian. 
(c) The bias dependence of the cut-off frequency for the low frequency Lorentzian. (d) The gate 
dependence of the cut-off frequency for the low frequency Lorentzian. Blue points correspond to 
the negative slope 8I/dQ g > (see fig. [2]). Red points correspond to the positive slope dI/dQ g < 
(see fig. [2]). The error bars are extracted from the fits to the Lorentzians. The continuous lines 
(blue, red) show the bias and gate dependencies for the Lorentzian cut-off frequency calculated in 
the described model. 
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positive slope (dI/dQ g > 0) (red curve in figure HJ), the situation is different. For negative 
bias voltage the Lorentzian cut-off frequency is continuously growing from f co ~ 60 kHz and 
reaches a maximum (f co ~ 120 kHz) close to zero bias voltage. For the positive bias voltage 
it rapidly decrease from the maximum to the initial value f co ~ 60 kHz). Figure Hfb), shows 
the gate dependence for the Lorentzian cut-off frequencies for both slopes (dI/dQ g <0). As 
is clearly shown in this figure, the gate dependence for the positive and negative slopes are 
different. The gate dependence for the positive slope (dI/dQ g > 0) (red curve) has a peak, 
with a small negative offset on the gate charge from the SET open state. The Lorentzian 
cut-off frequency behavior on the negative slope (8I/dQ g < 0) (blue curve) is a step like 
function. 

By integrating the Lorentzians in the noise spectra (see figure [3(b)) we can calculate the 
total variation of induced charge on the SET island for both fluctuators. The variation of 
induced charge, for the low-frequency fluctuator is of the order of Aqif ~ 6.6me rms . The 
same estimation for the high-frequency fluctuator gives Aq^j rs 30me rms . 

V. DISCUSSION 

In this section we analyze the bias and the gate dependence of the cut-off frequency of 
the high frequency Lorentzian (f co ~ 50 — 150kHz). In particular, we try to explain why 
the cut-off frequency is different for different biasing conditions. 

In our analysis we have assumed that there are in principle two possible sources for 
the low frequency noise, resistance fluctuators or charge fluctuators. The physical nature 
of the resistance fluctuators is not well understood, but they can be related with charge 
fluctuations. For instance, a charge oscillating in the tunnel barrier may modify both the 
transparency and the induced island charge. 

The noise from a resistance fluctuator in one of the tunnel barriers would have an asym- 
metry along the onset of the SET open state, as it was shown in section [TT] (see fig. H(b)). 
In order to explain the bias dependence of the cut-off frequency (see fig. IH(a)) in terms 
of resistance fluctuators we must assume that there is an individual resistance fluctuator 
located in each of the SET tunnel barriers. Furthermore these fluctuators must have the 
same tunneling rates. With this strong assumption, however, it is impossible to explain the 
sharp drop of the experimentally measured gate dependence of the cut-off frequency (see. 
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fig- SKb)). 

Thus, in order to explain the results for the high-frequency Lorentzian, we will assume 
that there are individual charge ffuctuators affecting the SET, and that each Lorentzian in 
the experimentally measured spectra is due to a single ffuctuator coupled to the SET island. 
Many experimental groups have suggested a microscopic nature of these ffuctuators. The 
microscopic nature is not known well, but there are suggestions that it could be traps in the 
substrate dielectric close to the SET or in the aluminum surface oxide. 

Here we will argue that the sources of these two level ffuctuators are located outside the 
barrier and that they may have a mesoscopic nature. In reference^ it is argued based on 
electrostatic analysis of the tunnel barrier, that such ffuctuators could not be localized inside 
barrier. There are also other experiments^ 1 ^, where it is argued that the charge ffuctuators 
most probably are localized outside the tunnel barrier. 

A typical SET is made from thin aluminum films which are not uniform; they consist of 
small grains connected to each other. In figure [5^a) we are show a SEM image of a sister 
sample to the measured one. In figure [5](b) we also show an AFM image of the same sister 
sample. It can be clearly seen that there are many small grains close to the device. We will 
assume that some grains are separated from the main film by a thin oxide layer but also 
capacitively connected to the SET island. 

Electrostatically such a grain can be described as a single electron box—, capacitively 
coupled to the SET island with capacitance C g and having a tunnel contact with resistance 
R' and capacitance C with one of the bias leads as indicated in figure[5]^c). The situation is 
almost equivalent if the grain would be tunnel connected to the SET island and capacitively 
connected to the SET bias lead. For a detailed analysis we should estimate the energy scales 
for this grain. We assume that the linear dimension of the stray aluminum grain is in the 
range 1 — 5 nm. The charging energy for this grain is of the order E' c = e 2 / (2C^) ~ 10 _1 eV, 
where = C + C' g + C' env , and C' env is the capacitance to the rest of the environment. This 
charging energy is substantially larger than the experimental temperature and the charging 
energy of the SET (ksT <C Eq <C E' c ). In addition there will be be further separation of 
the energy levels due to the small size of the grains. 

The ratio of capacitances C' g /C^ is given directly by the charge induced on the SET 
island which we already have extracted by integrating the Lorentzians. Thus for the high 
frequency Lorentzian we have C' g /C^ = 0.030 and for the low frequency Lorentzian this ratio 
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FIG. 5: (a) A SEM image of a sister sample to the measured one. The black bar shows 100 nm 
linear scale, (b) An AFM pictures of the edge of an aluminum film on the SiOx surface (c) 
Equivalent electrostatic scheme, where the small metallic grain is capacitively coupled to the SET 
island and has a tunnel connection with a bias lead. 



is 0.0066. 

In our model the single electron box is capacitively coupled to the SET island, and tunnel 
coupled with one of the SET leads. The average potential of the SET island acts, in this 
system, as a gate potential for the single electron box and induces the charge n' g = C' g 4>/e 
on the grain. The charging dynamics for the electron box can be described by the orthodox 
theory using a master equation approach^. Electron tunneling changes the number of excess 
electrons n in the grain. The differences in the electrostatic energy, when electrons tunnel 
to (+) and from (— ) the grain are: 



A£±(n) = 2E' C (±n ^n' g + l/2). 



(2) 



The tunnel rates of electron tunneling to or from the grain is a function of the tunnel 
resistance R' and the electrostatic energy gain AS^(n): 

1 AS±(n) 



w; 



e 2 R' 1 - exp (-A£±(n)/(AkT)) ' 
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(3) 



The probability a n to have n excess electrons in the grain obeys the master equation: 
da r . 



(H K-i a n-i + w n+1 a n+1 - a n (w+ + w n ) (4) 

For our case of low temperature, ksT <C E' c there are only two nonvanishing probabilities, 
a n and o n +i. This simplifies the problem, and it is convenient to define the distance from the 
grain degeneracy point, i.e. the point where these two nonvanishing states are degenerate, 
as An' g = C' g ((j) — 4>o)/e. Here 0o = e(n + 1/2) jC' g is the SET island potential needed to 
reach the grain degeneracy point. 

In the time domain, the dynamics for charging this grain from the lead is a random 
telegraph process, and the spectrum of this process is a Lorentzian function with a cut-off 
frequency defined by the sum of charging and escaping rates^: 

lS.nl f F' \ 

U = w+ + w~ +1 = -^L coth I ^An;J = A{<f> - O ) coth {B{<j> - O )) , (5) 

where A = C' g /{eC^R') and B = eC g /(2k B TC' s ). 

From eq. ([5]) we thus see that the cut-off frequency of the Lorentzian depends directly 
on An' g and therefore on the potential of the SET island, (p. When the grain is biased away 
from its degeneracy point, i.e. when An' g > 2k B T '/ 'E' c , the cut-off frequency grows linearly 
with the SET island potential. This means that if we are far from the grain degeneracy 
point, the cut-off frequency is relatively high and the relative frequency shift due to the 
change in will be small. On the other hand if we are close to the grain degeneracy point 
the cut-off frequency will be close to its minimum and the relative change in frequency due 
to can be substantial. The maximum relative frequency change occurs when the potential 
just barely reaches the grain degeneracy point and can be calculated from eq.0 Taking in 
to account the bounds of the island potential, — e/(2Cs) < <t> < e/(2Cs), we get 

fco,mi 

We note that the relative frequency shift is independent of R' and that a large charging 
energy of the SET is important to observe a frequency shift. Thus it is clear that the relative 
frequency shift of most Lorentzians will be very small. To observe a frequency shift as large 
as in figureHl the the grain will have to be close to its degeneracy point, and in addition 
the charging energy of the SET will have to be large, to create a substantial frequency 
shift. Obviously a relatively strong coupling between the grain and the SET island is also 
important. 
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FIG. 6: (color online) (a) The SET island potential calculated from the orthodox theory, as a 
function of bias voltage and gate charge, (b) Line cuts along the two diagonals in (a), where the 
measurements have been performed. As can be seen the potential is asymmetric with respect to 
the SET degeneracy point. If the grain SEB is biased away from its degeneracy point the frequency 
of the tunneling on and off the grain should be proportional to the island potential. The resulting 
fit is shown in figure[71 

The average potential of the SET island will depend both on the bias and gate voltages 
of the SET, 4>(V, V g ), and can be calculated using the orthodox theory^. We have calculated 
the potential as a function of bias voltage and gate voltage, as is shown in figureEta), and 
in particular we have calculated the potential along the two diagonals (roughly where the 
measurements have been performed) which is shown in figure[6]^b) to give an idea of how the 
potential varies (c.f. figuresd] and EJ^a)). As can be seen, the potential is asymmetric with 
respect to the SET degeneracy point and varies between — e/(2Cs) and e/(2Cs). 
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FIG. 7: (color online) The cut off frequency of the high frequency Lorentzian versus the potential 
for the SET island. The potential is calculated from the orthodox theory for the bias and gate 
voltage at which the spectrum was recorded. Stars (red) represent measurements taken at positive 
gain (D points), whereas circles (blue) represent measurements taken at negative slopes (B points). 
The line is a fit to eq.[5j 

Inserting the calculated potential into eq. ((Sj) we can thus make a quantitative comparison 
between our model and the measured data. In figured we plot the cut-off frequency of the 
high frequency Lorentzian as a function of the SET island potential. As can be seen there 
is a good agreement between the data and our model, (eq.EJ). To obtain this we have used 
three fitting parameters, namely A, B and <fto of equation (J5j). From this fit we can extract 
important parameters of the grain SEB. From the A parameter we get the tunnel resistance 
to the grain, using the capacitance ratio extracted from the integrated charge change of 
the Lorentzian Aq^f. The extracted tunnel resistance is R' = 2.4 GQ. Considering the 
size of the grain and that the sample was exposed to air before the measurement, this is 
quite reasonable. From the B parameter we can extract the temperature of the electrons in 
the lead, T\ ea ^ = 130 mK. The experiment is performed at a mixing chamber temperature 
of 25 mK, however we also have to take into account that SETs are always substantially 
overheated above the cryostat base temperature^. Typically the overheating of an SET 
island is a few hundred mK. The lead next to the island will be colder, which is consistent 
with the temperature we extract. We can also see from the fit that the potential of the SET 
island passes through O = 0.22, where the grain is at its degeneracy point. 

From the above discussion we can now also consider the low frequency Lorentzian and try 
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to understand why its cut-off frequency is not changing with bias or gate voltage. As can be 
see in figureH^c and d) we do not observe any significant change with the island potential. 
Assuming that the low frequency Lorentzian is also due to a grain, then this fluctuator must 
be far from its grain degeneracy point and therefore the relative frequency change is very 
small. In this case the tunnel resistance to this grain must be substantially higher than that 
of the grain responsible for the high frequency Lorentzian. 

If we consider a large ensemble of many grains close to the tunnel junctions of an SET, 
most of them will only be weakly coupled to the island and will thus together make up a 
1/ / background in the noise spectrum. A few of the grains may be more strongly coupled 
to the island and will show up as individual Lorentzians in the noise spectrum just as we 
and many others have observed. In general it will be relatively rare that a strongly coupled 
grain is close to its degeneracy point since the charging energy of the grain is much larger 
than that of the SET. In addition we do need a large charging energy of the SET in order 
to change the frequency of the fluctuator. Thus it will be quite rare that a Lorentzian will 
show the large change in the cut-off frequency which we have observed. Therefore it is not 
so surprising that we have not been able to find this type of behavior in more samples in 
spite of numerous efforts. 

The results presented here does not exclude that there are also other types of fluctuators 
with different physical mechanisms that contribute to the noise, however we are confident 
that this is one type of fluctuator which we have been able to identify and found a good way 
to characterize. 

VI. CONCLUSION 

In conclusion we have measured the noise of a single electron transistor from a few Hz up 
to 10 MHz. We find a spectrum with two Lorentzians superimposed on a 1/f background. 
The cut-off frequency of one of the Lorentzians depends strongly on the bias and gate voltages 
whereas the other does not. Our data is consistent with a model where the low-frequency 
noise comes from the random charge process of two effective single electron boxes coupled 
to the SET. We suggest that these single electron boxes are due to small aluminum grains 
coupled by tunneling to one of the leads and capacitively to the SET island. We are able to 
fit our data to this model with good agreement and we can extract parameters for the one 
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of the fluctuators. 
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